Mitochondrial genomes revisited: why do different lineages retain different genes?

The mitochondria contain their own genome derived from an alphaproteobacterial endosymbiont. From thousands of protein-coding genes originally encoded by their ancestor, only between 1 and about 70 are encoded on extant mitochondrial genomes (mitogenomes). Thanks to a dramatically increasing number of sequenced and annotated mitogenomes a coherent picture of why some genes were lost, or relocated to the nucleus, is emerging. In this review, we describe the characteristics of mitochondria-to-nucleus gene transfer and the resulting varied content of mitogenomes across eukaryotes. We introduce a 'burst-upon-drift' model to best explain nuclear-mitochondrial population genetics with flares of transfer due to genetic drift.

A DRP-like pseudoenzyme coordinates with MICOS to promote cristae architecture.

Mitochondrial cristae architecture is crucial for optimal respiratory function of the organelle. Cristae shape is maintained in part by the mitochondrial inner membrane-localized MICOS complex. While MICOS is required for normal cristae morphology, the precise mechanistic role of each of the seven human MICOS subunits, and how the complex coordinates with other cristae shaping factors, has not been fully determined. Here, we examine the MICOS complex in Schizosaccharomyces pombe, a minimal model whose genome only encodes for four core subunits. Using an unbiased proteomics approach, we identify a poorly characterized inner mitochondrial membrane protein that interacts with MICOS and is required to maintain cristae morphology, which we name Mmc1. We demonstrate that Mmc1 works in concert with MICOS complexes to promote normal mitochondrial morphology and respiratory function. Bioinformatic analyses reveal that Mmc1 is a distant relative of the Dynamin-Related Protein (DRP) family of GTPases, which are well established to shape and remodel membranes. We find that, like DRPs, Mmc1 self-associates and forms high molecular weight assemblies. Interestingly, however, Mmc1 is a pseudoenzyme that lacks key residues required for GTP binding and hydrolysis, suggesting it does not dynamically remodel membranes. These data are consistent with a model in which Mmc1 stabilizes cristae architecture by acting as a scaffold to support cristae ultrastructure on the matrix side of the inner membrane. Our study reveals a new class of proteins that evolved early in fungal phylogeny and is required for the maintenance of cristae architecture. This highlights the possibility that functionally analogous proteins work with MICOS to establish cristae morphology in metazoans.

Evolutionary History of Oxysterol-Binding Proteins Reveals Complex History of Duplication and Loss in Animals and Fungi.

Cells maintain the specific lipid composition of distinct organelles by vesicular transport as well as non-vesicular lipid trafficking via lipid transport proteins. Oxysterol-binding proteins (OSBPs) are a family of lipid transport proteins that transfer lipids at various membrane contact sites (MCSs). OSBPs have been extensively investigated in human and yeast cells where 12 have been identified in Homo sapiens and 7 in Saccharomyces cerevisiae. The evolutionary relationship between these well-characterized OSBPs is still unclear. By reconstructing phylogenies of eukaryote OSBPs, we show that the ancestral Saccharomycotina had four OSBPs, the ancestral fungus had five OSBPs, and the ancestral animal had six OSBPs, whereas the shared ancestor of animals and fungi as well as the ancestral eukaryote had only three OSBPs. Our analyses identified three undescribed ancient OSBP orthologues, one fungal OSBP (Osh8) lost in the lineage leading to yeast, one animal OSBP (ORP12) lost in the lineage leading to vertebrates, and one eukaryotic OSBP (OshEu) lost in both the animal and fungal lineages.

The persistent homology of mitochondrial ATP synthases.

Relatively little is known about ATP synthase structure in protists, and the investigated ones exhibit divergent structures distinct from yeast or animals. To clarify the subunit composition of ATP synthases across all eukaryotic lineages, we used homology detection techniques and molecular modeling tools to identify an ancestral set of 17 ATP synthase subunits. Most eukaryotes possess an ATP synthase comparable to those of animals and fungi, while some have undergone drastic divergence (e.g., ciliates, myzozoans, euglenozoans). Additionally, a ∼1 billion-year-old gene fusion between ATP synthase stator subunits was identified as a synapomorphy of the SAR (Stramenopila, Alveolata, Rhizaria) supergroup (stramenopile, alveolate, rhizaria). Our comparative approach highlights the persistence of ancestral subunits even amidst major structural changes. We conclude by urging that more ATP synthase structures (e.g., from jakobids, heteroloboseans, stramenopiles, rhizarians) are needed to provide a complete picture of the evolution of the structural diversity of this ancient and essential complex.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria).

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Evolution: Divergent trajectories predate the origins of animals and fungi.

Animals, fungi, and their closest protist relatives comprise the clade Opisthokonta. Although they are comparatively closely related, animals and fungi have diverged greatly from one another. A new study demonstrates that the genomic features that are characteristic of animals and fungi arose even before the origin of these two kingdoms.

Single cell genomics reveals plastid-lacking Picozoa are close relatives of red algae.

The endosymbiotic origin of plastids from cyanobacteria gave eukaryotes photosynthetic capabilities and launched the diversification of countless forms of algae. These primary plastids are found in members of the eukaryotic supergroup Archaeplastida. All known archaeplastids still retain some form of primary plastids, which are widely assumed to have a single origin. Here, we use single-cell genomics from natural samples combined with phylogenomics to infer the evolutionary origin of the phylum Picozoa, a globally distributed but seemingly rare group of marine microbial heterotrophic eukaryotes. Strikingly, the analysis of 43 single-cell genomes shows that Picozoa belong to Archaeplastida, specifically related to red algae and the phagotrophic rhodelphids. These picozoan genomes support the hypothesis that Picozoa lack a plastid, and further reveal no evidence of an early cryptic endosymbiosis with cyanobacteria. These findings change our understanding of plastid evolution as they either represent the first complete plastid loss in a free-living taxon, or indicate that red algae and rhodelphids obtained their plastids independently of other archaeplastids.

First report of mitochondrial COI in foraminifera and implications for DNA barcoding.

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in many marine environments and play a major role in global carbon cycling. Species recognition in Foraminifera is mainly based on morphological characters and nuclear 18S ribosomal RNA barcoding. The 18S rRNA contains variable sequence regions that allow for the identification of most foraminiferal species. Still, some species show limited variability, while others contain high levels of intragenomic polymorphisms, thereby complicating species identification. The use of additional, easily obtainable molecular markers other than 18S rRNA will enable more detailed investigation of evolutionary history, population genetics and speciation in Foraminifera. Here we present the first mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences ("barcodes") of Foraminifera. We applied shotgun sequencing to single foraminiferal specimens, assembled COI, and developed primers that allow amplification of COI in a wide range of foraminiferal species. We obtained COI sequences of 49 specimens from 17 species from the orders Rotaliida and Miliolida. Phylogenetic analysis showed that the COI tree is largely congruent with previously published 18S rRNA phylogenies. Furthermore, species delimitation with ASAP and ABGD algorithms showed that foraminiferal species can be identified based on COI barcodes.

Depletion of a Toxoplasma porin leads to defects in mitochondrial morphology and contacts with the endoplasmic reticulum.

The voltage-dependent anion channel (VDAC) is a ubiquitous channel in the outer membrane of the mitochondrion with multiple roles in protein, metabolite and small molecule transport. In mammalian cells, VDAC protein, as part of a larger complex including the inositol triphosphate receptor, has been shown to have a role in mediating contacts between the mitochondria and endoplasmic reticulum (ER). We identify VDAC of the pathogenic apicomplexan Toxoplasma gondii and demonstrate its importance for parasite growth. We show that VDAC is involved in protein import and metabolite transfer to mitochondria. Further, depletion of VDAC resulted in significant morphological changes in the mitochondrion and ER, suggesting a role in mediating contacts between these organelles in T. gondii. This article has an associated First Person interview with the first author of the paper.

A Eukaryote-Wide Perspective on the Diversity and Evolution of the ARF GTPase Protein Family.

The evolution of eukaryotic cellular complexity is interwoven with the extensive diversification of many protein families. One key family is the ARF GTPases that act in eukaryote-specific processes, including membrane traffic, tubulin assembly, actin dynamics, and cilia-related functions. Unfortunately, our understanding of the evolution of this family is limited. Sampling an extensive set of available genome and transcriptome sequences, we have assembled a data set of over 2,000 manually curated ARF family genes from 114 eukaryotic species, including many deeply diverged protist lineages, and carried out comprehensive molecular phylogenetic analyses. These reconstructed as many as 16 ARF family members present in the last eukaryotic common ancestor, nearly doubling the previously inferred ancient system complexity. Evidence for the wide occurrence and ancestral origin of Arf6, Arl13, and Arl16 is presented for the first time. Moreover, Arl17, Arl18, and SarB, newly described here, are absent from well-studied model organisms and as a result their function(s) remain unknown. Analyses of our data set revealed a previously unsuspected diversity of membrane association modes and domain architectures within the ARF family. We detail the step-wise expansion of the ARF family in the metazoan lineage, including discovery of several new animal-specific family members. Delving back to its earliest evolution in eukaryotes, the resolved relationship observed between the ARF family paralogs sets boundaries for scenarios of vesicle coat origins during eukaryogenesis. Altogether, our work fundamentally broadens the understanding of the diversity and evolution of a protein family underpinning the structural and functional complexity of the eukaryote cells.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans.

BACKGROUND: The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered. RESULTS: We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome. CONCLUSIONS: Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

A functional bacteria-derived restriction modification system in the mitochondrion of a heterotrophic protist.

The overarching trend in mitochondrial genome evolution is functional streamlining coupled with gene loss. Therefore, gene acquisition by mitochondria is considered to be exceedingly rare. Selfish elements in the form of self-splicing introns occur in many organellar genomes, but the wider diversity of selfish elements, and how they persist in the DNA of organelles, has not been explored. In the mitochondrial genome of a marine heterotrophic katablepharid protist, we identify a functional type II restriction modification (RM) system originating from a horizontal gene transfer (HGT) event involving bacteria related to flavobacteria. This RM system consists of an HpaII-like endonuclease and a cognate cytosine methyltransferase (CM). We demonstrate that these proteins are functional by heterologous expression in both bacterial and eukaryotic cells. These results suggest that a mitochondrion-encoded RM system can function as a toxin-antitoxin selfish element, and that such elements could be co-opted by eukaryotic genomes to drive biased organellar inheritance.

Constructive Neutral Evolution 20 Years Later.

Evolution has led to a great diversity that ranges from elegant simplicity to ornate complexity. Many complex features are often assumed to be more functional or adaptive than their simpler alternatives. However, in 1999, Arlin Stolzfus published a paper in the Journal of Molecular Evolution that outlined a framework in which complexity can arise through a series of non-adaptive steps. He called this framework Constructive Neutral Evolution (CNE). Despite its two-decade-old roots, many evolutionary biologists still appear to be unaware of this explanatory framework for the origins of complexity. In this perspective piece, we explain the theory of CNE and how it changes the order of events in narratives that describe the evolution of complexity. We also provide an extensive list of cellular features that may have become more complex through CNE. We end by discussing strategies to determine whether complexity arose through neutral or adaptive processes.

Independent accretion of TIM22 complex subunits in the animal and fungal lineages.

Background: The mitochondrial protein import complexes arose early in eukaryogenesis. Most of the components of the protein import pathways predate the last eukaryotic common ancestor. For example, the carrier-insertase TIM22 complex comprises the widely conserved Tim22 channel core. However, the auxiliary components of fungal and animal TIM22 complexes are exceptions to this ancient conservation. Methods: Using comparative genomics and phylogenetic approaches, we identified precisely when each TIM22 accretion occurred. Results: In animals, we demonstrate that Tim29 and Tim10b arose early in the holozoan lineage. Tim29 predates the metazoan lineage being present in the animal sister lineages, choanoflagellate and filastereans, whereas the erroneously named Tim10b arose from a duplication of Tim9 at the base of metazoans. In fungi, we show that Tim54 has representatives present in every holomycotan lineage including microsporidians and fonticulids, whereas Tim18 and Tim12 appeared much later in fungal evolution. Specifically, Tim18 and Tim12 arose from duplications of Sdh3 and Tim10, respectively, early in the Saccharomycotina. Surprisingly, we show that Tim54 is distantly related to AGK suggesting that AGK and Tim54 are extremely divergent orthologues and the origin of AGK/Tim54 interaction with Tim22 predates the divergence of animals and fungi. Conclusions: We argue that the evolutionary history of the TIM22 complex is best understood as the neutral structural divergence of an otherwise strongly functionally conserved protein complex. This view suggests that many of the differences in structure/subunit composition of multi-protein complexes are non-adaptive. Instead, most of the phylogenetic variation of functionally conserved molecular machines, which have been under stable selective pressures for vast phylogenetic spans, such as the TIM22 complex, is most likely the outcome of the interplay of random genetic drift and mutation pressure.

Homologue replacement in the import motor of the mitochondrial inner membrane of trypanosomes.

Many mitochondrial proteins contain N-terminal presequences that direct them to the organelle. The main driving force for their translocation across the inner membrane is provided by the presequence translocase-associated motor (PAM) which contains the J-protein Pam18. Here, we show that in the PAM of Trypanosoma brucei the function of Pam18 has been replaced by the non-orthologous euglenozoan-specific J-protein TbPam27. TbPam27 is specifically required for the import of mitochondrial presequence-containing but not for carrier proteins. Similar to yeast Pam18, TbPam27 requires an intact J-domain to function. Surprisingly, T. brucei still contains a bona fide Pam18 orthologue that, while essential for normal growth, is not involved in protein import. Thus, during evolution of kinetoplastids, Pam18 has been replaced by TbPam27. We propose that this replacement is linked to the transition from two ancestral and functionally distinct TIM complexes, found in most eukaryotes, to the single bifunctional TIM complex present in trypanosomes.

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists.

Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.

Neutral evolution of cellular phenotypes.

Eukaryotes exhibit a great diversity of cellular and subcellular morphologies, but their basic underlying architecture is fairly constant. All have a nucleus, Golgi, cytoskeleton, plasma membrane, vesicles, ribosomes, and all known lineages but one have mitochondrion-related organelles. Moreover, most eukaryotes undergo processes such as mitosis, meiosis, DNA recombination, and often perform feats such as phagocytosis, and amoeboid and flagellar movement. With all of these commonalities, it is obvious that eukaryotes evolved from a common ancestor, but it is not obvious how eukaryotes came to have their diverse structural phenotypes. Are these phenotypes adaptations to particular niches, their evolution dominated by positive natural selection? Or is eukaryotic cellular diversity substantially the product of neutral evolutionary processes, with adaptation either illusory or a secondary consequence? In this paper, we outline how a hierarchical view of phenotype can be used to articulate a neutral theory of phenotypic evolution, involving processes such as gene loss, gene replacement by homologues or analogues, gene duplication followed by subfunctionalization, and constructive neutral evolution. We suggest that neutral iterations of these processes followed by entrenchment of their products can explain much of the diversity of cellular, developmental, and biochemical phenotypes of unicellular eukaryotes and should be explored in addition to adaptive explanations.

ER-shaping atlastin proteins act as central hubs to promote flavivirus replication and virion assembly.

Flaviviruses, including dengue virus and Zika virus, extensively remodel the cellular endomembrane network to generate replication organelles that promote viral genome replication and virus production. However, it remains unclear how these membranes and associated cellular proteins act during the virus cycle. Here, we show that atlastins (ATLs), a subset of ER resident proteins involved in neurodegenerative diseases, have dichotomous effects on flaviviruses-with ATL2 depletion leading to replication organelle defects, and ATL3 depletion to changes in virus production pathways. We characterized non-conserved functional domains in ATL paralogues and show that the ATL interactome is profoundly reprogrammed following dengue virus infection. Screen analysis confirmed non-redundant ATL functions and identified a specific role for ATL3, and its interactor ARF4, in vesicle trafficking and virion maturation. Our data identify ATLs as central hubs targeted by flaviviruses to establish their replication organelle and to achieve efficient virion maturation and secretion.

A single-cell genome reveals diplonemid-like ancestry of kinetoplastid mitochondrial gene structure.

Euglenozoa comprises euglenids, kinetoplastids, and diplonemids, with each group exhibiting different and highly unusual mitochondrial genome organizations. Although they are sister groups, kinetoplastids and diplonemids have very distinct mitochondrial genome architectures, requiring widespread insertion/deletion RNA editing and extensive trans-splicing, respectively, in order to generate functional transcripts. The evolutionary history by which these differing processes arose remains unclear. Using single-cell genomics, followed by small sub unit ribosomal DNA and multigene phylogenies, we identified an isolated marine cell that branches on phylogenetic trees as a sister to known kinetoplastids. Analysis of single-cell amplified genomic material identified multiple mitochondrial genome contigs. These revealed a gene architecture resembling that of diplonemid mitochondria, with small fragments of genes encoded out of order and or on different contigs, indicating that these genes require extensive trans-splicing. Conversely, no requirement for kinetoplastid-like insertion/deletion RNA-editing was detected. Additionally, while we identified some proteins so far only found in kinetoplastids, we could not unequivocally identify mitochondrial RNA editing proteins. These data invite the hypothesis that extensive genome fragmentation and trans-splicing were the ancestral states for the kinetoplastid-diplonemid clade but were lost during the kinetoplastid radiation. This study demonstrates that single-cell approaches can successfully retrieve lineages that represent important new branches on the tree of life, and thus can illuminate major evolutionary and functional transitions in eukaryotes. This article is part of a discussion meeting issue 'Single cell ecology'.